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FINITE ELEMENT ANALYSIS OF A COMPOSITE 
MATERIAL INTERFACE 

By 

Kenneth H. Murray* 

SUMMARY 

A finite element model of a composite material interface is developed 
to study the influence of the interface on the thermal strain in the composite. 
A plane stress model is used with an axisymmetric model as a check. The 
interface thickness, thermal coefficient, modulus, Poisson's ratio and the 
percent of mineral in the composite are variables in the study. The 
results confirmed the usability of the finite element model in studying 
the polymer-mineral interface. 


INTRODUCTION 

Composite materials can be considered to be a combination of two separate 
phases or materials with an interphase or interface between them. The inter- 
face transmits the forces and displacements between the phases and thus 
makes possible the combined action of the phases[|FT]. These forces and 
displacements have been analyzed by means of finite elements to estimate 
the matrix stresses [2-4] as well as the effect of the interface [5-7]. 

In all of these analyses, perfect bonding and zero thickness interface 
between the phases was assumed, although such an assumption is physically 
unreal istic. 

Consequently, a finite element analysis was made of a composite material 
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consisting of a continuous polymer phase {the matrixi) containing a discon- 
tinuous mineral phase (the filler). The phases were separated by a finite 
thickness interface to which a range of properties could be assigned. The 
effects of the interface properties on the thermal strain of the composite 
were determined and are reported in the next section. 

MODELS AND ANALYSIS 
Composite Model 

The model polymer-m.ineral composite material had properties which 
approximated those of an epoxy-glass system in which the gdass is in the 
form of small, uniformly sized, beads. The properties were: 


Polymer Mineral 

Elastic modulus, GN/m^ 3.047 68.95 

Thermal expansion coef. , K“^ 88.6x10"^ 8.5x10"^ 

Poisson's ratio 0.44 0.22 

Specific gravity 1.15 2.45 

Bead diameter, fm — 125 


A composite Wttth these properties was being investigated concurrently in an 
experimental program at NASA Langley Research Center. 

Finite Element Model 

The finite element model of the composite was basically a plane stress 
model of a unit cell of polymer surrounding a single glass bead, separated 
by a finite thickness interface (Figure 1). The bead was modeled as a 
hexagon because only straight line boundaries were permitted with the plane 
stress elements available in the computer. program which was used in the 
analysis. The interface, which can vary in thickness and material properties 
is essentially a geometrical extension of the hexagon. The bead diameter 
was taken as the distance between opposite sides of the hexagon. Since 
the diameter was fixed (125 ym), different volume percentages of mimeral 
were obtained by changing the size of the unit cell. Thus a large mineral 



volume percentage (e.g. 60%) required a smaller unit cell than did a small 
percentage (e.g. 5%). 

Using the basic unit cell, four finite element models were constructed. 

The first model (Figure 2) had a flexible link element as the interface. 

This link has been used successfully to model the bonding between concrete 
and reinforcing steel [8]. The material properties of this interface element 
were independent of the polymer and mineral because of the linking arrangement, 
so these properties could be varied over an extremely wide range. The second 
model used the same geometry except for the interface element (Figure 3). 

The interface element was a normal element geometrically but had its material 
properties bounded by the material properties of the polymer and mineral; 
e.g. it is unreasonable to expect the interface modulus to be larger than 
the mineral or smaller than the polymer modulus. The thi»id model was made 
up of two of the second models stacked on top of each other (Figure 4). 

This was done to determine lif the strain for two unit cells would differ 
from that for one cell. The fourth model was an axisymmetric idealization 
of the bead and polymer (Figure 5). This model used the same type of interface 
element as the second model. The idealization showhnin Figure 5 was rotated 
360° about the z-axis to develop a three-dimensional solid. 

Analysis 

The four models were analyzed using the finite element method by way 
of NASTRAN. NASTRAN, an acronym for MSA Stru ctural ^alysis, is a general 
purpose finite element computer program for structural analysis. A discussion 
of the finite element method and of NASTRAN is given in Appendix A. 

The loading system for the models was a temperature change from 300 K 
(approximately room temperature) to 400 K. In the finite element method 
the thermal displacements are found by defining a set of nodal loads which 
approximate the thermal load on an element. The nodal loads from each 
element are summed to obtain the total load vector. This approximation 
along with the basic idealization of polymer-mineral composite are necessary 
to analyze the system using finite elements. 

The boundary conditions for the first three models were the same 



(Figures 3 and 4). The x-axis nodes (16, 17, 21, and 22) were restrained 
in the y-direction and the y-axis nodes (18, 22, and 23) were restrained 
in the x-direction. The boundary conditions approximated those of a pdiece 
of the composite two unit cells wide, resting on a plate, with the composite 
axis of symmetry co/imcident with the y-axis. As a result of the symmetry 
only half of the piece (one unit cell) had to be analyzed. The boundary 
conditions for the fourth (axi symmetric) model were zero displacement along 
the r-axis in the y-direction (see Figure 5). 

The variables used in the analysis are listed in Table I. One composite 
variable, the volume percent of mineral, was included to confirm the results 
obtained with the other variables; if the mineral volume is increased, the 
thermal strain should decrease [9]. Four interface variables were used, 
the interface thickness plus the interface properties of elastic modulus E, 
coefficient of thermal expansion a, and Poisson's ratio p. The interface 
properties were assumed to change littdie over the 300-400 K temperature 
range investigated, and indeed, the large drop in the epoxy modulus does 
not occur until a temperature of 440 K is reached [10]. The interface prop- 
erties were assigned values equal to the polymer, the mineral, or the average 
of the two values. The particular value is indicated in the text and the 
tables by a lower case suffix p (polymer), m (mineral), or a (average). 
Thus the designation aa, Em means that the average value of the coefficiient 
of thermal expansion (48. 55xlO“®/K) and the mineral value of the elastic 
modulus (68.95 GN/m?) were assigned to the interface. Only the interface 
properties were varied; the polymer and mineral always were assigned their 
respective values listed in Table I. 

The geometry of the models, the variable mineral volume, and the 
variable interface properties created a major data preparation problem. 

To reduce the effort and the possible errors associated with the input data, 
a computer program was written to punch the NASTRAN input deck. This program 
is given in Appendix B along with an outline of its operation. 

RESULTS AND DISCUSSION 


The results of the finite element program (NASTRAN) are displacements 



at the nodes and stresses within the element. The displacements were divided 
by the length (height or width) of the unit cell to obtain strain values. 

These strain values are listed in Tables II through VII. An analysis of 
variance (ANOVA) of selected data is given in Tables VIII through X. 

Preliminary Results 

Model One: No results are presented in this report because the model 

never gave reasonable results. The material properties of the linking element 

20 

were varied systematical ly from zero to 10 and no combination of values 
gave reasonable results. 

Poisson's ratio: Table II shows some of the results when the Poisson's 

ratio was varied in Model Two. It is obvious from the table that this 
quantity does not affect the thermal strain values, therefore it was eliminated 
from further study and an average value was used (y = 0.33) for the interface. 

Interface thickness: Most researchers agree that the interface is a 

thin layer between the mineral and the polymer. Using this, the larger 
value of F did not appear valid. The remainder of the study uses only 
F = 1 and 10 ym. 

Model Three: The two cell model gave strain values identical to the 

one cell model (Model Two). The mineral content was 34% and the interface 
was 10 ym with all combinations of am, ap, Ep, and Em. 

Variables examined: The final results are based on the following 

variables with their discrete values listed in parentheses. All combinations 
were studied for: 

Volume Percent of Mineral - PER (5, 34, and 60%) 

Interace thickness - F (1 and 10 ym) 

Interface thermal coefficient - a (ap, am, and aa) 

Interface modulus - E (Ep, Em, and Ea) 

Model sstwo and four 

The reference model was Model Two with 34% mineral and an interface 
thickness of 1 ym. 
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Final Results 


Mineral content: Comparison of data in Tables III, IV, and V and the 

<( 

ANOVA in Table VIII indicatestitbattitbemtfiitineralcGonteiiitiiisaavyepyssignif (leant 
variable. The thermal strain in the model decreases with increasing mineral 
content. This can be predicted because the mineral is less responsive to 
temperature change than the polymer, =thus the more mineral the less thermal 
strain. 


Interface thickness: Comparison of datac.in Tables IV and VI and the 


ANOVA i n Tabl e IXtaridi Gatesli.thatfel;heiiiiintefefaGfel«'thi ckness 3i s^^a^sigh f iccint- vari abl 


Model: Comparison of data in Tables IV and VII and the ANOVA in Table X 

for Models Two and Four indicates no apparent significance for the model. 

This, of course, is good because now only the plane stress model needs to be 
studied in detail. 

Strain location: Comparison of Tables III through VII and the ANOVA 

in Tables VIII and IX indicates that the strain location is a significant 
variable. This is somewhat disturbing because it should not be significant. 

The variationsinhstran'hriailionggthbebbundciKyyis plotted in Figure 6. The low 
strain values in the middle region near the mineral bead are caused by the 
mineral's presence. 

Table XI gives the results of the basic model with the boundary conditions 
varied at the corner node (22). Standard conditions fix the node in both the 
X- and y-di recti ons. To help explain the honJuniformsstraninaacrosstthe-tOp 
and side of the model, this node (22) was released first in the x- and then 
in the y-di recti on. Only the two outside nodes and node 22 are presented, 
to observe the difference between them. It is quite obvious that the boundary 
conditions affect the final strains. The only way to eliminate this effect 
is to increase the number of elements near the boundaries to isolate the effect. 

Interface material properties: Observation of Table IV for Model Two, 

F = 1 ym and 34% mineral, indicates that as a and E increase the thermal 
strain increases. The thermal coefficient a has the most effect as can 
be seen from ANOVA Table 'IXII.; ItTsho'ucTdTibebinoitedethatathfehtotalaef’fect . 
of both variables is less than 1% in general. This presents critical 
questions as to the real effect of the interface on the final thermal strains. 



CONCLUDING REMARKS 


The basic procedure used in this study was to model the polymer- 
mineral composite as a plane stress elasticity problem modeled with finite 
elements. A unit cell with one bead and a variable interface was idealized 
as plane stress triangles and quadrilaterals. The interface was varied 
in thickness and material properties through a realistic range to obtain 
data on its effect on the thermal strain of the unit cell during a temperature 
change. 

For the standard model used in this study (34% mineral, interface thickness 
approximately il.% of the bead diameter;i, 1 ym, and interface material properties 
equal to the polymer, mineral and their average) the following conclusions 
can be drawn: 

1. The major factors influencing the thermal strain of the composite are: 

a) Mineral content 

b) Interface thickness 

c) Interface thermal coefficient 

2. The thermal strain decreases with increasing mineral content. 

3. The thermal strain decreases with increasing interface stiffness. 

4. More elements and nodes are needed to eliminate the variation in 
displacement if the plane stress elements are used. 

5. Experimentaili results are needed to establish the correct interface 
properties. 

6. The plane stress finite element idealization does a reasonable job 
'ipepreSenttingiithetpolymer.T.miTnefailireomposiinte'jtbutbneedsesomerrimore, e 
refinement. 

7. A three-dimensional representation needs to be made to evaluate the 
usefulness of the plane stress model beyond what Was done here. 
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Appendix A 


Finite Element Analysis - NASTRAN 

I. Basic Procedure (Reference A.l) 

A general continuous body can be idealized using the 
finite element method by dividing it into a discrete number 
of small elements. These elements are interconnected only 
at a finite number of points called nodes or joints. At each 
node a finite number of degrees of freedom, usually displace- 
ments and rotations, are prescribed. Thus the real system 
with displacements continuous throughout the body is idealized 
by a discrete system with displacements specified only at node 
This approach to structural analysis has been used on framed 
structures for quite some time, but in the last twenty years 
it has been extended to include two- and three-dimensional 
continua. Because it is a systematic procedure it is very 
well suited for implementation on digital computers. It is 
this fact plus the versatility of the method that has led to 
its great popularity. 

The first step in the analysis procedure is to idealize 
the continuum as a discrete number of finite elements. For 
each element a functional relationship between the internal 
displacements and the nodal displacements is written. Next 
the stress-strain-displacement relationships are developed 
and used to write the potential energy of the element. From 
the potential energy equation the stiffness matrix for the ele 
ment is derived. The stiffness matrices of all elements are 



iA-'2 

then summed to determine the equations of equilibrim for the 
continu\am. These equilibrium equations are then solved for 
the nodal displacements using modern numerical techniques. 

Once the nodal displacements have been computed, the forces-, 
stresses and strains in each element are found. 

II. Thermal Stress Analysis 

The temperature load is input as a nodal force which is 
equivalent to the thermal strain associated with the tempera- 
ture change. The thermal stress analysis by the finite element 
method can be best summarized by the following steps: 

1) Idealization 

The continuum is subdivided into a system of discrete or 
finite elements, interconnected at a discrete number of nodes. 
The unknown parameters of the problem are normally the dis- 
placements or rotations at each node. 

2) Functional 

A displacement function is selected which expresses the 
internal displacement of the element in terms of the unknown 
nodal displacements: 

{f} = [N] {6}g Cl) 

3) Strain 

The state of strain within an element is expressed in 
terms of the nodal displacements : 

{£} = [BJ {6}g . C2) 
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4) Stress 

The state of stress within an element is expressed in 
terms of the nodal displacements using the element strains, 
initial strains and constitutive properties of the material. 

In equation form: 

{a} = [D] ({e} - {Eq}) , (3a) 

or 

{a} = [D] ([B] {6}g - {£q}) . C3b) 

5) Total Potential Energy 

The total potential energy is expressed in terms of the 
nodal displacements using the stress and strain functions: 

n = U + V (4) 

where U is the internal strain energy 

U = 1/2 {e}'^ {a} dV_ 

e 

or after substituting equations (2) and (3b) 

U = 1/2 {6}^ [B]'^[D] [B]{6}gdVg-l/2/^ {6}^ [B] ^ [D] { }dV^ . (5) 

e ® 

The potential energy of the applied forces can be written as 
V= /^- {f}'^{p} dVg - /g{f}'^{g} dS - (6) 
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where {p} = body force per unit volume 

{g} = surface tractions per unit area 
S = surface (area) over which {g} acts 
{Fj^} = applied nodal forces 

For the case of thermal strain V = 0 and is omitted from the 
following. The total potential energy then is U,; . 

n = 1/2 /g {6 }^[B]T[d] [B]{6}^dVg - 1/2/ {6}T[B]T[D]{e^}dVe (7) 

6) Minimize Total Potential Energy 

The equilibrium of the system can be established by 
minimizing the total potential energy of the system. In equa- 
tion form this can be stated: 



e 


Performing the minimization, after using equation (1), yields 

/ [B]^[D] [B] {6} dV^ - [B]'^[D]{e_}dV = 0. (8) 

V e s Vq o e 

e ^ 

7) Nodal Forces 

In the minimization above four load types are usually 

- 

present, but here only one is present. Each must be converted 
into equivalent nodal forces. The second integral in equation 
(8) represents nodal loads resulting from initial strains or 

{Fq} = / [B]^[DJ {e^ldVg . (9) 

''e 
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The equivalent nodal force then is 

{F}^ = {F^} (10) 

8) Stiffness Matrix 

Another outgrowth of the minimization is the basic 
element relationship between load and displacement or the 
stiffness matrix [k], . s The basic equilibrium equation in 
discrete cooredinates can be written as 

{F}^ = IkU6}g 

and comparing it to equation (8) defines 

[kj = [B]'^[D] [B]dVg . (11) 

9) Assemblage 

Each element has now been reduced to a set of equilibrium 
equations. The continuum can now be represented by the sum of 
the equilibrium equations. If N denotes the total number of 
elements , then 

{F} = [K]{6} C12) 

N 

where {F} = Z {f} U3a) 

^ e 

e=l 

N 

[K] = E [k] (13b) 

e=l ® 

N 

and {6} = {6} . (13c) 

— T e 
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10) Boundary Conditions 

The external constraints are applied using modern numer- 
ical technique iprocedures . 

11) Solution 

The unknown displacements of equation (12) with con- 
straints are obtained using appropriate numerical procedures 
to solve the equations . 

12) Stress Recovery 

The displacements are used to determine the stress and 
strains in the continuum using equations (2) and (3a) . 

This entire procedure has been programed many times. The 
program used in this study is NASA's Structural Analysis Pro- 
gram - NASTRAN. The input for NASTRAN must contain the fol- 
lowing information; 

1) Type of analysis (static) 

2) Type of element (plane stress) 

3) Nodal point locations 

4) Element connections 

5) Material properties 

6) Loads (temperature change) 

7) Boundary conditions 

8 ) Output requirements 

Most of the input is quite easy, but the nodal locations 
and element connection cards can be difficult. For this reason 
a FORTRAN program was written to generate the geometry of the 
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polymer-mineral composite and then convert this into NASTRAN 
input format for punching. Appendix B gives a detailed de- 
scription of this procedure. 
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Appendix B 


NASTRAN Input Generation Program 

A FORTRAN program was written to generate the geometry 
of the mineral polymer composite and punch it in NASTRAN for- 
mat. The Executive Control and Case Control decks are read 
into this program as data for the final NASTRAN deck. Also 
read into this program are the GRDSET, PROPERTY, MATERIAL and 
TEMP cards from the Bulk Data deck. All other cards are 
generated and punched. The mineral bead diameter, mineral 
content in percent volume , interface thickness and the number 
of basic blocks in the model are also read into this program. 
This allows great flexibility in defining the variables of this 
polymer-mineral composite study. 

The program is listed on the following pages with a 
typical data deck. 
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DIMENSION deck- (30c 18 ) *GR I DX ( 35 ) *GR I DY ( 35 ) • IE ( 22 J ♦ JE ( 22 ) »XE ( 2 2) 
DIMENSION LE(22) 

B * BEAD RADIUS IN MICRONS 

F = INTERFACE THICKNESS 

PER * PERCENT BY VOLUME OF BEADS 

A 3 INSCRi*BED CIRCLE OF LARGE HEX 

R e RADIUS TO node OF LARGE HEX 

IN = LOGICAL UNIT FOR INPUT 

100 = LOGICAL UNIT FOR PUNCHED OUTPUT 

IC = LOGICAL UNIT FOR PRINTED OUTPUT 

NO » NUMBER OF HEXES 

MAXIMUM OF 2 FOR J^ODEL ONE AND 1 FOR MODEL TWO 
NDEC = NUMBER OF LEADING CARDS 
BHEX = RADIUS OF BEAD HEX 
ITYPE = TYPE OF INTERFACE 

0 FOR MODEL ONE AND 1 FOR MODEL TWO 

SET UP INPUT-OUTPUT UNITS 

IN*2 

XO-3 

100=7 

READ (IN*1J B.F*PER»NO*NDEC tITYPE 
1 FORMAT (3F10«0»312) 

NDECX=NDEC+1 . 

READ IN NASTRAN EXECUTIVE CONTROL DECK* 

CASE CONTROL DECK AND 
THE FOLLOWING BULK DATA CARDS 
GRDSET 
PROPERTY 
MATERIAL 
TEMP# 

WRITE U0#600) 

600 FORMAT (IHl) 

DO 15 I«l#NDECl 

RFAD (IN.n)(DECK(I*J)*J = l*18l 
WRITE! 10*11) (DECK( I *J) #J=1»18» 

15 CONTINUE 
11 FORMAT (18AA) 

BHEX = 1*15469*B 
R * SQRT(75.0/PER)*B/.86603 
CRIDX(l) * 0* 

GRIDY(l) * 0* 

IF (ITYPE) 300*19*300 

LINKAGE ELEMENT MODEL ONE 

19 DO 20 1=2*13 

20 GRIDX(I)«BHEX 
DO 21 1=14*25 

21 GRIDXT I)»BHEX+F 
DO 22 1=26*31 

22 GRIDX(I)-R 

JJ-1 


1-3 

GRinY(2)*0 

23 GRIDY(I)«GRIDY(I-1)4-60* 

IF 24#30»30 

2«^ IF t JJ-5 ) 25%lfi%n 

25 I-l+l 
JJ-JJ+1 
GO TO 23 

26 I-I+l 

JJ-1 

GRIDY(n-0. 

I-t+l 

GO TO 23 

27 WRITE(IO*3» 

3 FORMAT ( *Pf<OBLEMS») 

STOP 

30 DO AO I*»l*3l • 

RR-GRIDX( I ) 

THa(-GRIOY( I) +60.) *0*01745329 
GRIDXn )=RR*COStTH)+R 
40 GRIOY( I )=RR*S1N(TH)+0.86603*R 
CRIDX(32)-2.*R 
GR1DX(33)*’.*R 
GRIDX{34)»0* 

GRIDX(35»*0. 

GRIDY{32)-1*73206*R 
GFIDY<33)«0* 

GR10Y(34)«0* 

GRIDY(35)-1*73206*R 
DO 50 I-l»6 
IE(I)*l 
JE(n«3+l-l 
LEdUO 

50 KE(I)»2+I-1 
JE(6)=2 
DO 60 Is7*12 
IEm*2+I-7 
JE(I)»9+I-7 
KE(n-21*«-I-7 
60 LE(I)» 14*1-7 
JE(12)>8 
KE(12J*20 
DO 70 1=13*18 
IE(I )«20*I-13 
JE(n=21*I-13 
XEn)«27*I-13 
70 LE(IU26+I-13 
JE(18)»20 
XE(18)«26 
DO 80 1=19*22 
IE{I)»27*I-19 
JEd)=32+I-19 
XE(I)=26*I-19 
80 LE(n«0 
Jn(21)*30 
JE(22»*31 
<F(21)*34 
XE(22)»35 
DO 16 I=1,NDEC 
16 WRITE( 100*11) (OECK( I »J) *J=1*18) 
CD 90 1=1*35 



n n n 


90 WRITE (100*4) I*GRIDX(I)#GRIDYm 

4 FORMAT ( 'GRID* *4X*I8»8X#2F8*3I 
316 DO 100 1*1*6 

100 WRITE (100*5) I*IE(I)*JE(I)*KE(n 

5 FORMAT { •CTRMEM»2X*!e*7X**l»*3I8) 

DO no I»19*22 

no WRITE ( 100 * 6 ) I*lE(!)*JE(I)*KE(n 

6 FORMAT ( 'CTRMEM* *2X*I8»7X* *2' #318) 

DO 120 1=7*12 

120 WRITE (100*7) I. IE(I)*JE(I)*KE(I)*LE(II 

7 FORMAT ( »CODMEM' *2X*I8*7X* *2» *418) 

DO 130 1=13*18 

130 WRITE (100*8) I * I E ( I ) * JE ( I ) *XE ( I) *LE ( I) 

8 FORMAT ( ‘CQDMEM* *2X* I8*7X* *1* *418) 

IF ( 1 rYPE)380*139*380 

SPC CARDS 

139 11*1 
12*34 
13*12 
14*30 
I5»l 

WRITE! 100*12) 11*12*13*14*15 
12 FORMAT! •SPC*5X*3I8*8X*2ia» 

11*2 

12*35 

13*1 

14*29 

15*2 

WRITE! 100*12)11*12*13*14*15 

11*3 

12*28 

13*2 

14*33 

WF.ITE( 100*12)11*12*13*14*15 

149 IF(NO-1)27*240*140 

140 DC 141 1=1*35 

141 GRIDY( I )=GRIDY( I )+l*73206»R 
DO 150 1*1*27 

I1*I^35 • 

150 WRITE! 100*4) II *GRIDX( I ) *GRIDY(I) 

11*63 

WR1TE(IOO*4)II*6RIDX(30) *GRIOY(30) 
ll*6^> 

WRITE! I00*4)1I*GRIDX(31) *GRI0Y(31) 

11*65 

WRITE! 100*4) 11 *GR1DX( 32) *GRIOY( 32) 
11 * 66 - 

WRITE! 100*4)11*GRIDX(35) *GR1DY( 35) 

DO 151 1*1*6 
IE!I)*36 
JE( I )«38-*I-l 
XE!I)*37-fI-l 

151 LE(I)*0 • 

JE(6)*37 

DO 16') 1*7*12 
lE(I)-37+I-7 
JE!I )-44*-I-7 
KE!n»56-H-7 



u u u 


160 Lr.(n«49^i-7 
JE(l2»-<^3 
<F<12J«55 
DO 170 I«13*18 
IE< n=*56+I-l3 

170 LE(H»55-H-13 
I£a 8)=55 
J£(13)=62 
JE(14)=26 
JEn5)»31 
JEn6)»63 
JE(17)»64 
JEa8l»6l 
KE(13»=6l 
KE(14)^62 
KE(15)«26 

ICE(16)=31 • 

KEa ?)-63 
KE(ir)=64 
DO IbO I»19«22 

180 LE(I)«0 
IE<19)*62 

IE(20)»26 , 

IEUD-31 
IE<22)*63 
JE<19»«65 
JE(20)»32 
JE(21)*63 
JE(22)*64 
XE <19 1*61 
KE(20M62 
KE(2l»a35 
ICE (22) *66 
DO 200 I«l*6 
11*1*22 

200 WRITEn00t5)ntIE(I)»JE(I)tlCEm 
DO 210 1*19.22 
Il*I+22 

210 WRITEn00.6)Il.IE(I).JE(I)«ICE(I) 

DO 220 1*7.12 
Il*I-^22 

220 WRITE(I00.7)1I.IE(I) .JE(I) .ICE(I)>LE(I) 

00 230 1*13.18 
11*1+22 

230 WRlTE<I00.8)II.IE(I).JE(n«ICE(n»U£(X) 

11*4 

12*63 

13*1 

U.66 

15*1 

SPC CAROS 

WRITE( 100.12) II « 12. 13.14.15 
IF<NO-l)27.240.24l 

240 WRITE( 100.13) 

13 FORMAK 'SPCADD’7X.*100».7X.*1*.7X.*2* .7X.»3*) 
GO Tt 131 

241 WRiTtC 100.14) 

14 FORMAT! *SPCADD»7X* 100 ‘7X*1*7X*2 *7X* 3 *7X*4» ) 

131 WRITE (100.9) 


ri n n non 


9 FORMAT (‘ENDDATA*) 

STOP 

REGl^LAR ELEMENT MODEL TWO 

300 00 310 I=2»7 

310 GRlDXtn=BHEX 
DO 311 I=8tl3 

311 6RIDX( DaBHEX-t-F 

DO 312 Ial4*l9 

312 GRIDX(I)aR 

JJ=l 

I»3 

6RIDY(2)*0 

323 GRIDYn)=GRIDYn-l)460« 

IF (1-19)324*330*330 

324 IF (JJ-5) 325*326*327 

325 I = I<*1 • 

JJaJJ+l 

GO TO 323 

326 I=I+l 

JJ»1 

GR|DY(I)*0« 

I*I + l 
GO TO 323 

327 WRITE(IO*3l 
STOP 

DEVELOPS AND PUNCH GRID AND CONNECTION CARDS 

330 CD340 1*1 *19 
RRaGRIOX(I) 

TH=(-GRI0Y(I )+60.)*0. 01745329 
GFIDX(I)=RR*COS(TH)-*-R 
340 GRIDY( I )aRR*SIN(TH)+0#86603*R 
GFIDX(20)=2«*R 
GRIDX(21)*2»*R 
GRIDX(22)»0« 

GRIDX(23)=0. 

GRIDY(20)=1*73206*R 

GRIDY(21)=0« 

GRIDY(22)=0. 

GRIDV(23)=1«73206*R 
D0350 I»l*6 
IEm«l 
JE(n-34l-l 
LE(I)-0 

350 KE(I)a2+I-l 
JE(6)a2 
00 360 I»7#18 
IE(I)»2-^I-7 
JErn»3+I-7 
KE(I)a94l-7 
360 LE{I)a8+l-7 
JEC12»«2 
JE(18)-8 
XEn2)*8 
XE(18)*14 
DO 370 2-19*22 
KE(n-20+I-19 
370 LE(n=0* 


IE(19)-14 
IE(20)*15 
IE(2X)-17 
Ict22)*18 
JE(19)«15 * 

JE(20)*16 
JE(21)-18 
JE(22)-19 
D0416 I=1»N0EC 

416 WRITE( lOOtll) (0ECIC(I*JUJ«1*18» 

D0390 1=1*23 

390 WRITE (100*4) I *GR IDX ( I) »GRIDY( 1 ) 

GO TO 316 
C 

C DEVELOPE AND PUNCH SPC CARDS 

C 

380 11*1 
12*22 
13*12 
14*18 
15*1 

WRITE( 100*12) 11*12*13*14*15 
11*2 . 

12*23 

13*1 

14*17 

15*2 

WRITE (100*12) 11 *12* 13 *14, 15 

11*3 

12*16 

13*2 

14*21 

WRITE( 100*12) 11*12*13*14*15 
GO TO 149 
END 
C 

C EXAMPLE DATA DECK FOLLOWS 

C 

/* ' 
125* 1*0 5* 126 1 

ID POLYMER*MINERAL 
•.PP DIE? 

SOL 1*0 
TIME 10 
CEND 

TITLE » ONE HEX - 5 PERCENT MINERAL 

subtitle = INTERFACE = 1 * VARYING BOND REGULAR ELEMENT 

FCHO = BOTH 

bPCFORCES * ALL 

OLOAD = ALL 

ELSTRESS * ALL 

DISPLACEMENT * ALL 

SUBCASE 1 

LABEL = ALfHA * MIN * E = MIN AND MU * 0*33 
TEMPERATURE ( LOAD) *1 
SPC *100 
BEGIN BULK 

GRDSET 3456 

PTRMEM I I !• 

PTRMEM 2 2 1« 


22 


PODMEM 

1 2 

!• 


p::dmem 

2 3 

1« 


MATl 

1 68*95 

• 22 

1 

MATt 

2 3*047 

*44 

1 

MATl 

2 68*95 

• 33 

X 

TEMPO 

1 100* 




8o5E-6 

88*6E-6 

8.5E-6 


TABLE I 

Variables and their Values 


Mineral Content, 
by Volume 

PER 

5, 

34 and 60% 

Interface Elastic 
Modulus 

E 

Ep 

Ea 

Em 

(3.047 GN/m^) , 
(35.9985) and 
(68.95) 

Interface Thermal 
Coefficient 

a 

ap 

aa 

am 

(88.6 X 10"^/°K) 
(48.55 X 10“® and 
(8.5 X 10-6) 

Interface Thickness 

F 

1, 

10, and 20ym 

Interface Poisson's 
Ratio 

y 

yp 

and 

(0.44), ya (0.33) 
ym (0.22) 

Unit Cell Length 
X-Direction 

hx 

(1118.0, 428.7, 322.7)ym 

Unit Cell Length 
Y-Direction 

Ly 

(968.2, 371.3, 279.5)ym. 
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TABLE II 




Effect of Poisson's Ratio 




Average % 

Strain 

For F = lOym 




ym 


ya 

yp 

am 

Em 

0.5136 


0.5135 

0.5134 

ap 

Em 

6: 5603 


0.5621 

0.5603 

am 

Ep 

0.5249 


, 0.5244 

0.5234 

ap 

Ep 

0.5518 


0.5538 

0.5556 



Average % 

Strain 

For F = lym 




ym 


ya 

yp 

am 

Em 

0.5491 


0.5491 

0.5491 

ap 

Em 

0.5537 


0.5540 

0.5543 

am 

Ep 

0.5502 


0.5502 

0.5501 

ap 

Ep 

0.5528 


0.5531 

0.5532 


Notes : 






All -values are % strain averaged over all 
exterior nodes 

Mineral content is 34% in model two 
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Node 22 y displacement 
Node 22 fixed 




Figure 1 Basic Unit Cell Model of Polymer r- Mineral Composite 



Polymer 



Figure 2 Model One 




Element 
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Figure 3 Model 


Figure 4 Model Three 



Figure 5 Model Four (Axisymmetric) 




